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We report scanning tunneling microscopy (STM) and spectroscopy (STS) measurements of mono- 
layer and bilayer WSe 2 . We measure a band gap of 2.21 =b 0.08 eV in monolayer WSe 2 , which is 
much larger than the energy of the photoluminescence peak indicating a large excitonic binding 
energy. We additionally observe significant electronic scattering arising from atomic-scale defects. 

Using Fourier transform STS (FT-STS), we map the energy versus momentum dispersion relations 
for monolayer and bilayer WSe 2 - Further, by tracking allowed and forbidden scattering channels as 
a function of energy we infer the spin texture of both the conduction and valence bands. We observe 
a large spin-splitting of the valence band due to strong spin-orbit coupling, and additionally observe 
spin-valley-layer coupling in the conduction band of bilayer WSe 2 - 
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Transition metal dichalcogenides (TMDs) have gained 
great popularity recently for their potential applications 
in optical and electronic devices pQ. The band gaps of 
these materials he in the visible part of the spectrum, 
and undergo an indirect-to-direct transition at monolayer 
thickness. In monolayers, the band extrema (valleys) lie 
at the K and K’ points [2], and the lack of inversion sym¬ 
metry leaves these valleys physically distinguishable (for 
example, via optical selection rules). The heavy tran¬ 
sition metals in TMDs also contribute to strong spin- 
orbit coupling (SOC) for band edge carriers, with espe¬ 
cially large effects expected for the tungsten dichalco¬ 
genides such as WSe 2 . This SOC leads to an out-of-plane 
Zeeman-type spin splitting of the valence band as large 
as 500 meV, with opposite spin polarizations at the K 
and K’ valleys due to time reversal symmetry leading to 
a coupling of the spin and valley degrees of freedom. 

AB-stacked bilayer TMDs are inversion symmetric, 
and as a consequence the spin and valley polarizations 
of the top layer are opposite those of the bottom, thus 
effectively restoring spin and valley degeneracy. They ad¬ 
ditionally exhibit an extra degree of freedom character¬ 
ized by a layer pseudospin. In WSe 2 , the SOC is larger 
than the interlayer hopping amplitude, thus localizing 
the carriers by layer near the Brillouin zone (BZ) cor¬ 
ners and leading to a coupling of the spin, valley, and 
layer degrees of freedom (see Ref. |T] for a review of these 
properties). Numerous spintronic and valleytronic appli¬ 
cations have been proposed to exploit these degrees of 
freedom in TMD monolayers and bilayers HJEH33]. 

Significant work has been performed to characterize 
these polarizations in WSe 2 with direct optical excita¬ 
tion [TOl [TT| 14. 15] and spin-resolved ARPES [16, ITT] . 
However, the former is unable to capture physics away 
from direct transitions (i.e. away from the BZ cor¬ 
ners) and the latter cannot probe unoccupied (conduc¬ 
tion band) states. Here we directly probe the electronic 
states of monolayer and bilayer WSe 2 via STM and STS 
measurements in order to gain a deeper understanding 
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FIG. 1. (a) Schematic of the experimental setup showing the 
STM tip and an optical microscope image of the measured 
WSe 2 on graphite sample, (b) Typical dl/dV spectroscopy 
of monolayer and bilayer WSe 2 on a log scale, (c) STM to¬ 
pography of a typical region of bilayer WSe 2 - There are two 
distinct defect types, with the large defects far less prevalent 
than the small defects. The scale bar is 20 nm. Inset: Atomic 
resolution of a small defect. The scale bar is 2 nm. 


of its band structure and internal quantum degrees of 
freedom. Fig. [lja) shows a schematic of our measure¬ 
ment setup. All measurements were performed in ultra- 
high vacuum at a temperature of 4.5 K. WSe 2 sits on a 
graphite flake to provide a conducting substrate for col¬ 
lecting the tunnel current. The sample was fabricated 
using a dry-peel transfer technique mem. WSe 2 and 
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graphite were mechanically exfoliated onto separate Si 02 
substrates, then the WSe 2 flake was picked up using a 
polycarbonate film and micro-mechanically aligned onto 
a graphite flake. Cr/Au electrodes were written using 
electron beam lithography, and the sample was annealed 
overnight in vacuum at 300 °C. 

Fig. 0b) shows representative dl/dV spectroscopy for 
monolayer and bilayer WSe 2 - We extract the electronic 
band gap of both by measuring the energy range where 
the dl/dV spectroscopy is zero, from which we find 
= 2.21 4= 0.08 eV for monolayer and 2.14 d= 0.05 eV 
for bilayer WSe 2 - This is significantly larger than the 
typically observed optical gaps in monolayer and bilayer 
WSe 2 (^1.65 eV and 1.6 eV, respectively) [2(71 [22] . in¬ 
dicative of a giant excitonic binding energy [23]. We 
note that this spectroscopically extracted gap actually 
represents the separation of the onset of tunneling into 
the valence and conduction band extrema, which, due to 
the finite width of the bands, is smaller than the sepa¬ 
ration between the center of the valence band maximum 
and the conduction band minimum. The peaks and/or 
inflection points in the dl/dV spectrum likely represent 
the band centers, which would imply a band gap a few 
hundred meV larger than that extracted above, in rea¬ 
sonably good agreement with the GW level DFT calcu¬ 
lated gap [24, [25]. Due to small changes in the tip work 
function resulting from the tip shaping procedure (for ex¬ 
ample, pulsing a gold particle onto the end of the tip), 
we have occasionally observed small (<150 meV) rigid 
shifts of the entire dl/dV spectrum in sample voltage. 
To account for this, we take dl/dV spectroscopy before 
all measurements, and when necessary offset the sample 
voltage to match the dl/dV reference spectra of Fig. 0b). 

Topographic maps of WSe 2 exhibit a significant den¬ 
sity of defects. Fig. 0c) shows one such example for bi¬ 
layer WSe 2 , in which two species of defects are present. 
There are small atomic-scale defects throughout the im¬ 
age (the inset is an atomically-resolved image of one such 
defect), along with two larger scale defects. These de¬ 
fects are more clearly visible in spatially resolved dl/dV 
maps, as shown in Figs. [2](a) and (b) for the same area 
as in Fig. 0c), taken at sample voltages of V s = -1.35 V 
and -1.5 V, respectively. In bilayer WSe 2 , defects of a 
given species appear with two clearly distinct strengths, 
where the weaker defects presumably reside in the lower 
WSe 2 layer. This indicates we are capable of probing 
electrons in both layers. In monolayer WSe 2 , all defects 
are observed with similar strength as they all reside in 
the same layer [26] . By counting the number of defects 
in numerous dl/dV maps taken on monolayer WSe 2 , we 
find a defect density of 1.1 db 0.3 x 10 12 cm -2 for the 
smaller defects, and approximately two orders of magni¬ 
tude lower density for the large defects. Given the similar 
appearance of sulfur vacancies observed in M 0 S 2 [27 it is 
likely that the small defects we observe here are due to se¬ 
lenium vacancies, although a full defect characterization 



FIG. 2. dl/dV maps at constant sample biases of (a) -1.35 V 
and (b) -1.5 V on bilayer WSe 2 taken over the same area 
as Fig. 0c). The two species of defects (large and small) 
are easily distinguishable. The small defects additionally ap¬ 
pear either strong or weak, depending on whether they are in 
the top or bottom layer. Friedel oscillations surround the de¬ 
fects, and their wavelength grows shorter as the energy moves 
deeper into the valence band. Scale bar is 20 nm for both. 
Insets: Fourier transforms of the maps. As the scattering 
wavelength grows shorter, the size of the scattering disk in 
momentum space grows larger. Scale bar is 3 nm -1 for both. 


is outside the scope of this work. 

These defects act as scattering centers for the elec¬ 
trons in WSe 2 , resulting in Friedel oscillations visible as 
concentric rings surrounding each defect [28 . The wave¬ 
length of these oscillations at a given energy is controlled 
by the band structure of the WSe 2 - So, for example, 
since the valence band grows wider as the energy is low¬ 
ered (Fig. 0a) to (b)), the wavelength of the oscillations 
grows shorter as they are inversely proportional. FT- 
STS provides direct visualization of the available scatter¬ 
ing channels for electrons in momentum space [29]. The 
long-wavelength Friedel oscillations of Figs. 0a) and (b) 
are due to intravalley scattering processes and show up as 
disk-like features at the center of their respective Fourier 
transforms (insets). The signal in the Fourier transform 
can be understood as a map of the joint density of states 
(JDOS) of the band structure. The intravalley scatter¬ 
ing process with the largest possible momentum transfer 
connects one side of a given valley to the other, setting 
the radius of the disk observed in the Fourier transform. 
The disk is then filled in by all other smaller momentum 
transfer scattering processes connected within a single 
valley. As a result, the radius of the central scattering 
disk gives a measure of the size of the Fermi surface at 
a given sample voltage. The blue dots in Figs. 0a) and 
(b) mark the momentum extracted from the size of these 
scattering disks as a function of sample voltage, providing 
experimental maps of the band structures of monolayer 
and bilayer WSe 2 (see the Supplemental Material [26] for 
details on wave vector extraction and band assignments). 
Data from different sets of dl/dV maps taken with dif¬ 
ferent tips and at different spots on the sample are all 
plotted together, as they follow the same general trend. 
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FIG. 3. Band structure of (a) monolayer and (b) bilayer 
WSe 2 - The solid red lines are ab initio calculations of the 
band structure. The calculated band gap is manually set to 
match the spectroscopy as described in the main text. Blue 
dots are the experimentally extracted momentum. Points are 
assigned to bands as described in the supplementary material. 
All points have an uncertainty in momentum of less than 10% 
of their pre-offset value (and typically less than 5%), arising 
from the uncertainty in identifying the size of the scattering 
disk in the Fourier transforms. There is an energy smearing 
between 5 and 20 mV on all points due to the ac voltage 
added to acquire the dl/dV signal. These have been omitted 
for clarity. 


Data points have been plotted on both sides of a given 
valley since the Fermi surfaces are roughly circular. 

We confirm our extracted energy versus momentum 
dispersion by comparing with the calculated ab initio 
band structures of monolayer and bilayer WSe 2 with 
spin-orbit coupling terms included (solid red lines in 
Fig. §, taken directly from Ref. [21J Since the size of 
the band gap is not well constrained in these calcula¬ 
tions, it has been manually set to match the band edges 
extracted from the dl/dV spectroscopy in Fig.[ljb) (plus 
a small offset to account for the finite band width). Our 
experimental results are in reasonable agreement with 
the theoretically anticipated bands; most notably the ef¬ 
fective masses generally match well with the theoretically 
predicted values. The extracted bands at T near -2 eV 
are around a hundred meV below the calculated bands 
for both monolayer and bilayer WSe 2 , suggesting possi¬ 
ble discrepancies between the DFT calculation and the 
true band structure. Disagreements may also result from 
mixing of the scattering signals of two bands. 

In addition to the intravalley scattering channels, there 
could, in principle, be scattering between different val¬ 
leys. Such scattering occurs with larger momentum 
transfer than intravalley scattering, and therefore the res¬ 
onances should be further from the center of the Fourier 
transforms. Due to strong spin-orbit coupling in WSe 2 , 
the valence band at the K point, which is primarily com¬ 
prised of m =±2 d-orbitals of W mm, is spin-split by 
around 450 meV, with bands having opposite spins in 
the K and K’ valleys. Fig. [4je) shows a schematic of 
the band structure of monolayer WSe 2 at an energy be¬ 


tween the two spin-split valence bands. There are six 
Fermi pockets surrounding the K and K’ valleys, alter¬ 
nating between spin up at K (green) and spin down at K’ 
(pink). Without this spin texture, atomic-scale defects 
would induce intervalley scattering, identified by reso¬ 
nances in the Fourier transform at 47r/3ao, where ao is 
the WSe 2 lattice constant (y/3 times closer to the cen¬ 
ter than the lattice resonances, and rotated from them 
by 30°). However, the spin polarization of the valleys is 
expected to suppress these resonances, as the intervalley 
scattering process now additionally requires a spin-flip 
which non-magnetic defects are unable to provide. 

We simulate the Fourier transform anticipated for 
each case by calculating the both the spin-dependent 
and -independent JDOS of the band structure shown in 
Fig. mm- The result is plotted in Fig.[4jf), where the 
colored features correspond to scattering resonances ex¬ 
pected only in the case of bands lacking spin polarization. 
The central resonance represents the intravalley scatter¬ 
ing discussed previously, and its diameter is twice the 
diameter of the K-point bands. Due to the geometry of 
the chords which can connect two points on the constant 
energy contour, the scattering signal decays away from 
the center of the Fourier transform but becomes large 
again just around the circumference of the disk, which is 
in good agreement with the features we observe in the 
maps shown in the insets of Fig. i Fig. gg) shows a 
comparable symmetrized Fourier transform of an atomi¬ 
cally resolved dl/dV map taken on monolayer WSe 2 [26) . 
The features anticipated from K-valley mixing are con¬ 
spicuously absent despite the high density of atomic-scale 
defects normally capable of inducing such scattering, im¬ 
plying the anticipated spin-polarization of the valleys. 
In contrast, this valley mixing is readily observable in 
graphene with atomic-scale defects [33], which lacks the 
spin polarization of the K-point valleys [34]. 

In the conduction band, the lowest energy bands ex¬ 
ist around the K and Q points, the latter of which is 
not a high symmetry point but occurs about halfway be¬ 
tween T and K. Fig. |4][a) shows a general schematic of 
the Fermi surface in the conduction band along with the 
spin polarization theoretically anticipated for WSe 2 m • 
While there is only a small spin-splitting anticipated at 
the K point (<50 meV) since these bands are predomi¬ 
nately m=0 d-orbitals of W [36] , the Q-point bands are 
expected to be split by over 200 meV in monolayer and 
500 meV in bilayer WSe 2 since they are comprised of 
both m —~tl p-orbitals of Se and m =±2 d-orbitals of 
W [301I3T] . Fig. ^b) shows the JDOS calculation for the 
conduction band with the same color coding scheme as in 
Fig- gf)- The K-point bands are removed from the cal¬ 
culation to reflect the fact that they are not expected to 
factor significantly into the observed scattering processes 
due to the reduced tunnel probability into K in the mono- 
layer [26] . and do not exist for the bottom-most ^250 
meV of the bilayer conduction band. Fig. [4][c) shows the 
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FIG. 4. (a) General schematic of the conduction band of monolayer and bilayer WSe 2 . The relative sizes of the bands depends 
on the layer number. In bilayer WSe 2 only the Q-point bands are present at the very bottom of the conduction band, and the 
K-point bands are spin-degenerate. The green denotes spin up bands and the pink denotes spin down, (b) Simulated JDOS 
of the band structure shown in (a). Only the Q-point bands are considered. Scattering resonances shown in white are always 
allowed, but colored features are only possible with no spin texture on the bands, (c) Experimental FT-STS on monolayer 
WSe 2 at V s = +1.5 V. The sharp resonances in a hexagon around the outside of the image are due to the atomic lattice. Only 
scattering allowed with the spin-textured bands is present (white features), implying the anticipated spin texture of the Q-point 
bands, (d) Experimental FT-STS on bilayer WSe 2 at V s = +1.3 V. The same subset of scattering features allowed in monolayer 
WSe 2 are observed, indicating strong layer polarization in bilayer WSe 2 . (e) Schematic of the band structure of monolayer 
WSe 2 between the spin-split valence bands, (f) Simulated JDOS of the band structure shown in (e). (g) Experimental FT-STS 
on monolayer WSe 2 at V s = -1.6 V. No intervalley scattering resonances are observed (blue features in (f)), implying the 
anticipated spin texture of the K-point bands, (h) Experimental FT-STS on bilayer WSe 2 at V s = -2.0 V. The intervalley 
scattering is also suppressed, similarly indicating strong layer polarization in bilayer WSe 2 - The scale bars are 10 nm -1 for all, 
and the white hexagons represent the first BZ. 


Fourier transform at the bottom of the conduction band 
of monolayer WSe 2 - Similar to the case of the valence 
band, only the features allowed assuming spin-polarized 
bands are present (note that the features corresponding 
to Q to Q scattering along the sides of the hexagonal 
Brillouin zone are present but very weak.) 

The spin degeneracy of the spin-split bands is restored 
in bilayer WSe 2 due to the inversion symmetry of the lay¬ 
ers. So, for example, if a given band in the bottom layer is 
spin down, that same band in the upper layer will be spin 
up. Despite the restored spin degeneracy, the Fourier 
transforms of our experimental dl/dV maps in the bi¬ 
layer conduction band (Fig. [4jd)) and the valence band 
(Fig. (4^h)) do not exhibit all the resonances anticipated 
in the simulation 133 • They do, however, match the exact 
subset of resonances expected given spin-polarized bands 
(white features in Figs.^b) and (e)). Put another way, 
the bilayer maps exhibit only the features observed in the 
monolayer maps, implying an effective spin-polarization 
of the bilayer WSe 2 bands. This spin-polarization arises 


from strong layer-polarization of the bilayer WSe 2 - Since 
there is only a small component of p z in the Q- and K- 
point bands, the spin-orbit coupling dominates over the 
interlayer hopping, thus enhancing the layer polarization. 
As we are capable of probing electrons in both layers of 
bilayer WSe 2 , the observation of nearly identical Fourier 
transforms for monolayer and bilayer WSe 2 provides di¬ 
rect evidence of strong spin-layer-valley entanglement in 
both the valence and conduction bands. 

We have presented local tunnel spectroscopy measure¬ 
ments of monolayer and bilayer WSe 2 - By examining 
electronic scattering from atomic-scale defects we are 
able to map out portions of the band structure of the ma¬ 
terial, as well as to identify scattering pathways forbidden 
by spin and layer polarizations of the bands. Our results 
directly demonstrate that intervalley electronic scatter¬ 
ing is suppressed in the valence band of monolayer WSe 2 
as a result of the large spin-splitting of those bands. Fur¬ 
thermore, we observe strong layer-polarization in bilayer 
WSe 2 , providing the first experimental evidence of the 
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spin-valley-layer coupling in the conduction band. While 
previous applications have been proposed to utilize this 
effect for slightly hole-doped bilayer WSe 2 [8], our re¬ 
sults present a path towards realizing these novel com¬ 
bined spintronic and valleytronic applications in devices 
of either doping. 
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